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Abstract. We consider an exclusion process representing a reactive dynamics 
of a pulled front on the integer lattice, describing the dynamics of first class 
X particles moving as a simple symmetric exclusion process, and static second 
class Y particles. When an X particle jumps to a site with a Y particle, their 
position is intechanged and the Y particle becomes an X one. Initially, there is 
an arbitrary configuration of X particles at sites . . . , —1, 0, and Y particles only 
at sites 1,2,..., with a product Bernoulli law of parameter p, < p < 1. We 
prove a law of large numbers and a central limit theorem for the front defined by 
the right-most visited site of the X particles at time t. These results corroborate 
Monte-Carlo simulations performed in a similar context. We also prove that the 
law of the X particles as seen from the front converges to a unique invariant 
measure. The proofs use regeneration times: we present a direct way to define 
them within this context. 



1. Introduction 

Few mathematical results exist about microscopic models of pulled front prop- 
agation representing non-equilibrium pattern formation in chemical reactions, or 
physical or biological phenomena (see [10] and [TT] for a review of the physical liter- 
ature). One-dimensional interacting particle systems which are microscopic versions 
of the Fisher-Kolmogorov-Petrovsky-Piscunov equation, representing evolutionary 
phenomena in genetics, have been studied in [T] and [3]. More recently, in several 
works (see [7J, [S] and [2]), systems of branching interacting random walks on the 
lattice were studied, which consider the spatial ordering of particles: two types of 
particles perform independent continuous time simple random walk movement: X 
particles which jump at rate Dx and Y particles at rate Dy. Upon contact with 
an X particle, a Y particle becomes X. In [3] the one-dimensional case where 
Dy = 0, representing the combustion of a propellant towards a stationary state, 
was analyzed. There, the particles are symmetric random walks with initially one Y 
particle at each site 1,2,... and an arbitrary configuration of X particles at —2,-1,0 
with a finite l\ norm with a certain exponential weight. If we call rt the position 
of the right-most visited site at time t by an X particle, it was proved that a.s. 
rt/t — > v, with v > and that e 1//2 (r e -i t — e~ 1 vt) converges to a Brownian motion 
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with non-degenerate variance. In [7], for the case of symmetric random walks with 
Dx = Dy a shape theorem was proved in arbitrary dimensions. In particular, in 
dimension d = 1 it was proved that if the initial configuration of all the particles is 
a product Poisson measure with a finite number of X particles, r% satisfies a strong 
law of large numbers. 

Mai, Sokolov, Kuzovkov and Blumen [8] performed Monte Carlo simulations for 
a variation of the above described model in which both the X and the Y particles 
perform symmetric simple exclusion. These numerical computations indicate that in 
the case Dx = Dy, with an initial condition which is a product Bernoulli measure, 
the front has a ballistic movement with normal fluctuations. 

In this paper we study a process where the X particles perform symmetric simple 
exclusion but where Dy = 0. We prove a law of large numbers and a functional 
central limit theorem for the position of the foremost visited site and for the number 
of activated particles, giving an indication that corroborates the behavior observed 
in the numerical simulations of [Sj. In the model we consider, there are two types 
of particles: the X particles which move as a symmetric simple exclusion process; 
the Y particles, which do not move. Initially there are no X particles at sites 
x > 0, while the configuration rj := {r/(0, x) : x < 0} of X particles at sites x < 0, 
where ry(0, x) is the number of X particles at site x, is such that 7/(0, 0) = 1 but 
otherwise arbitrary. Initially there are no Y particles at 0, —1, —2, . . . while at sites 
1,2,..., the Y particles are distributed according to a product Bernoulli distribution 
of parameter p. When an X particle jumps to a site where there is a Y particle, 
their position is interchanged and the Y particle becomes an X particle. Since the 
Y particles do not move, the dynamics can be defined in terms of the configuration 
of the X particles and the rightmost visited site at time t > 0, which we call rt- We 
adopt the convention that ro = 0. The state space of the process is then 

n := Ur,p,rj) : r G Z,p G Z,rj G {0, l}{--»- 1 ' r >J . 

Here p represents a counter for the number of Y particles which have been activated. 
The infinitesimal generator of the process is 

Lf(r,p,rj) := T, x , y < r ,\x-y\=iV(x)(^ ~ r)(y))(f{r,P,(Tx,yV) ~ f(r,P,v)) 
+pr)(r)(f(r + l,p+ l,ry + 5 r+1 ) - f(r,p,rj)) 
(1) +(1 - p)r](r)(f(r + l,p,rj) - f(r,p,rj)), 

where 5 X denotes the configuration with one particle at x, while (T xy r} denotes the 
configuration obtained from r\ after flipping the values of r\{x) and r\{y). 

We will use the notation rj(t, x) for the number of particles at time t > at site x, 
Pt the value of the counter for the number of activated Y particles (so that pt — po 
is the actual number of activated particles) and rt the position of the front at time 
t > 0. We will use the notation rj(t) := {rj(t,x) : x < rt} for the configuration 
of particles at time t and will call the process {(rt,pt,i](t)) : t > 0} the exclusion 
reactive process. 
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Theorem 1. Assume that initially r = p = 0, 7/(0,0) = 1 and n(0,x) € {0,1} is 
arbitrary for x < 0. 

(i) There exist av>0 and a w > which do not depend on the initial condition 
{n(0,x) : x < 0}, such that a.s. 

lim — = v, and lim — = w. 

t— >oo t t-^oo t 

(ii) There exist o\ and 02, < 01, 02 < 00, which do not depend on the initial 
condition {r/(0, x) : x < 0}, such that 



e x l 2 {r^ H - e~ x vt), and e 1/2 (p e ~i t - e~ l wt), t>0, 

converge in law as e —* to Brownian motions with variances o\ and 02, 
respectively. 

Theorem 2. Consider the process as seen from the front, r_ rt r/(i). There exist 
exactly two invariant measures: One supported on the configuration with no parti- 
cles, and another, /Joo . The domain of attraction of the first consists of exactly the 
configuration with no particles. Any nontrivial configuration in {0, l}^--^ 1 ' } is in 
the domain of the second; if we denote by fit the distribution of the process T^ n rj(t), 
then \it — * fJ-oo in the sense of weak convergence of probability measures. 

As in [2], this model does not satisfy any obvious sub-additivity property which 
would give a direct proof of part (i) of theorem[TJ But in a certain sense, the interac- 
tion given by the exclusion dynamics of this paper, is stronger than the independent 
random walk dynamics of [3], or the annihilating random walks of [2\. 

The proof of these results is based in regeneration time methods as discussed in 
[2] and [3]. Nevertheless, to reduce the tail estimates of these regenerations times to 
manageable expressions, we have mapped the exclusion reactive process to a zero- 
range reactive dynamics, with total jump rate 1 per site. This mapping and the 
definition of the regeneration times is given in section [2j 

In [2] and [3], the regeneration times where defined following [12], in terms of two 
alternating sequences of stopping times. One of these sequences defines the first 
times at which activated particles behind the leading one branch. This definition is 
performed using a space-time line which decouples the behavior of the old particles 
behind the front with respect to the behavior of the front itself. This is the approach 
used in this paper to define the regeneration times in section [2j However, we consider 
important to have an understanding from a more fundamental point of view about 
alternative ways which could be used to define regeneration times within the context 
of interacting particle systems representing reactive fronts. For this reason, in section 
El we present a definition of regeneration times for the exclusion reactive process, 
which is not done in terms of sequences of stopping times and which does not require 
a space-time line decoupling the dynamics at the left from the dynamics at the right 
of the front. The approach we present is in the spirit of Kesten [6] within the context 
of Random Walks in Random Environments. 
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In section El theorems [U and [2] are proved assuming that the regeneration times 
and the corresponding position of the front have finite second moments. In section 
[U these estimates are performed. The non-degeneracy of the variances o~\ and o~i of 
Theorem 1 is proved in section 

2. ZERO-RANGE REACTIVE PROCESS AND REGENERATION TIMES 

The first step in the proof of theorems Q] and [2j will be to couple the exclusion 
reactive process defined by ([I]), with a zero range process where particles branch 
at the right-most visited site. Then, it will be enough to prove a law of large 
numbers and a functional central limit theorem for the right-most visited site of this 
new process, together with the convergence towards an invariant measure for the 
law of this process as seen from its front. In the first subsection we will define the 
zero-range reactive process. In subsection l2.2| we will construct a version of the zero- 
range reactive process where particles are labeled. In subsection 12.31 we construct 
an auxiliary process, which gives lower bounds for the position of the front. Then, 
in subsection 12.41 we will define the regeneration times for the labeled version of the 
zero-range reactive process following some of the methods introduced in [6]. 

2.1. One dimensional zero-range reactive process. Consider a configuration 
(r,p,rj) of the stochastic combustion process with exclusion. Let q := r — p. We 
define £(<?) as the number of particles between the site r + 1 and the rightmost empty 
site. In other words, ((q) := r — x\, where x\ := sup{x < : r\[x) = 0}. Next define 
X2 := sup{x < x\ : r](x) = 0}, the position of the second right-most empty site, and 
C(q ~ 1) := x i ~ x 2 — 1) the number of particles between x% and x%. In general, for 
n > 2, we define x n := sup{x < x n _i : r](x) = 0}, while — n + 1) := x n -i —x n — \. 
Let := {C(x, t) : x < q} and qt : = r t —pt with initial condition (q^, (q) = (q, (). It 
is easy to check that the stochastic process {(qt,Pt, C(*))> * — 0} follows the dynamics 
of a one-dimensional reactive process with a zero-range dynamics with infinitesimal 
generator 

Lf(Q,P,0 = Z x ,y< q -.\ x -y\=lUv(x) >0)(f( q ,p,C x ' y )-f(q,P,0) 

+ (1 - p)l( V (q) > 0)(f(q + l,p,(-5 q + 5 q+l ) - f(q,p, ()) 

(2) +PUV(Q) > 0)(f(q,p + 1, C + S q ) - f(q,p, 0). 

This is a zero-range process with total jump rate 2 at those sites strictly to the 
left of the rightmost visited site q, with jump rate 1 to the left and 1 — p to the 
right at site q, and with branching at rate pl(C(q) > 0) at site q. Note that pt 
represents the number of times a branching has occurred. We will call the triple 
{(qt,Pt, C(0) • * — 0} the zero-range reactive process. 

2.2. Labeled process. We will make an explicit construction of the zero-range 
reactive process where each particle carries a label z£Z, representing the priority 
it has. The movement of a given particle will not be affected by particles with 
smaller labels. The construction will be performed in terms of a stochastic process 
{{y(f),qt) '■ t > 0}, where at time t > 0, the first component y(t) represents the 



FRONT PROPAGATION IN AN EXCLUSION DYNAMICS 



5 



positions on Z of a random number of particles, and qt the position of the rightmost 
visited site. Thus, the state space of this process is § := T x Z, where T := Uagz^ A - 
A typical element of T will be denoted by {y x : x £ A}, where A is the set of 
labels. We will furthermore use the notation p := sup{x : x £ A} and pt, when the 
corresponding set A t is time dependent. 

Let us fix an initial condition ({y x : x £ Ao},qo) £ §. Now associate to each 
x £ Aq L) {z £ Z : z > po} a discrete time simple symmetric random walk X x 
starting from 0, and a sequence {t x ^ : i > 1} of i.i.d. rate 2 exponential random 
variables, which will represent the potential jump times of an associated continuous 
time random walk. Let us also choose another sequence {vi : i > 1} of i.i.d. Bernoulli 
random variables of parameter p. We choose all these random variables independent 
of each other. 

Let us first define the dynamics of our process for an initial condition ({y x : x £ 
Ao},qo) such that Aq C Z is finite, and such that qo = sup{y x : x £ Aq}. We will 
associate to each discrete time random walk X x ,x £ Aq, a continuous time random 
walk Y x , such that Y x (0) = y x . Let n\ be the cardinality of the set {y x : x £ Aq}. We 
identify at each site in this set, the particle with the smallest label: let us call them 
xi, . . . , x ni . Let be the label within the group x±, . . . , x ni where the minimum 

7~i := min{rii\ ■ ■ ■ > } is achieved. If y Xt < qo, we let the continuous time random 
walk Y Xt jump at time t\ according to X Xt , while the other random walks do not 
move. Thus, the first change in the process {({Y x (t) : x £ A t },q t ) : t > 0} occurs at 
time 7~i when Y Xt (ri) = y Xt + X Xf (l). If y Xt = qo, and if the discrete time random 
walk X Xt jumps to the left at time 1, we let Y Xt jump to the left at time t±. Thus 
y x*{ti) = Vx* + ^x.(l)- If Ux, = Qo, but the random walk X Xlt jumps to the right 
at time 1, we let Y Xt jump to the right at time t\ only if the Bernoulli random 
variable v\ = 0, case in which Y Xm (r{) = y Xt + X Xt (l) and q Tl = qo + 1, while if 
v\ = 1, a random walk Y po+ i, which will follow the trajectory of the random walk 
X P()+ \ + qo, is created at time T\ while the remaining random walks do not move, so 
that A T1 = Aq U {po + 1}- This defines the dynamics of the process {(y(t),qt);t > 0} 
in the time interval [0, t±], with y(t) := {Y x (t) : x £ A t }.. 

Let us now recursively define the process for arbitrary times. Assume that for 
some k, such that k > 1, n& and have been defined and also the process in the 
time interval [0, r^]. Let pk ■= sup{x : x £ A Tk }. Call rik+i the cardinality of the 
set {Y x (rk) ■ x £ A Tk }, and identify at each site in this set the particle with smallest 
label: we call these labels x\, . . . ,x nk+1 . Let x* be the label where the minimum 

Tk+i := min{Tj.^ +1 ' ) , . . . , ri^ 1 ^} is achieved. Denote for 1 < i < as N Xi the total 
number of jumps performed up to time by the random walk Y Xi . If Y Xt (r^) < q Tk , 
we let the continuous time random walk Y Xt jump at time t^+i according to X Xt so 
that Y Xr (rk+i) = y x , + X Xt (N Xt + 1), while the other random walks do not move. 
If Y Xt (rk) = q Tk and if the random walk X Xt jumps to the left at the time N Xt + 1, 
we let Y Xt jump to the left at time 7>+i. If Y Xf (Tk) = q Tk but the random walk X Xt 
jumps to the right at the time N Xt + 1, we let Y Xt jump to the right at time t^+i 
only if the Bernoulli random variable «i+ Pfc - Po = 0, while if v i+ Pk - Po = 1 a random 
walk Y Pk+ i, following the discrete time trajectory of q Tk + X Pk+ i, is created. 
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Let us now consider the case in which Aq is not necessarily finite, but nevertheless 
po = sup{x : x G Ao} < oo. Consider the initial condition of positions of the particles 
y(0) = {i)x '■ x G Aq} and of the front qo. We use the following notations: given 
A C Z and n > 1, define y n (0) := {y x G y(0) : y x < q - n} and A n := {x G A : 
Vx < Qo~ n}. Since Aq is finite, we can define the process {(y n (t), qf) : t > 0} as in 
the previous paragraphs. 

Lemma 1. There exists a set of full measure such that for every t > the following 
statements are true. 

(i) There is an no such that if n > no, 

(ii) Let p" := sup{x : x G A™}. Then 

p t := lim p™, 

exists. 

(iii) For every x G At := ^4™°- 

Y x (t) := lim 

n— >oc 

exists. 

(iv) T/ie /imit 

q t := lim g t n , 

n— >oo 

exists. 

Proof. Let us prove part (i). Without loss of generality we assume that for each 
n > 1, Aq +1 ^ Aq. Consider the event E n := {A™ / A™ : for some m > n}. By the 
lemma of Borel-Cantelli it is enough to prove that 

oo 

(3) ^P[£ n ]<oo. 

71=1 

Now the probability of the event E n is upper bounded by the probability that in the 
time interval [0, t], a Poisson process of rate 2 has performed at least n steps. Indeed, 
the event E n is contained in the event that some particle initially at a distance larger 
than n from the front, is alone at the foremost visited site before time t. But this 
can happen only if at least n jumps where performed before time t. Thus, 

k=n 

which proves ([3]). Similar arguments can be used to prove parts (ii), (iii) and 
(iv). □ 

Define now y(t) := {(Y x (t) : x G At}. Then the triple {(y(t),q t ) : t > 0}, defines 
a probability measure P on the Skorokhod space D([0, oo); S), which we will call the 
labeled zero-range reactive process. It turns out that the particle count 
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C(t,x) = Yl UY x >(t)=x), 
x'eA t 

together with the pair p t and qt defined in lemma [TJ satisfies the dynamics defined 
by the infinitesimal generator ([2]). 

Consider an initial condition w := ({y x : x G Ao},qo) G §, such that y po = qo, 
where p = sup{x : x G A}. Define e) P0j(?0 G S as 5 pom := ({y x : x G B },q ), 
where Bq := {po}- This defines two coupled labeled zero-range reactive processes 
{{{Y x (t) : x G A t },q t ) : t > 0} and {({*£(*) : z G Sj.gQ : t > 0}, with initial 
conditions u; and <5 PO)(?0 respectively. The corresponding particle counts define two 
coupled zero range reactive processes with initial conditions corresponding to w 
and 5 POt q , which we will denote by ( w and Cs Po , qo respectively. The corresponding 

rightmost visited sites of these processes will be denoted by qf and q^ P0 ' q ° ; whereas 

the counters of the number of activated particles by pf and p t P °' q ° ■ Furthermore we 
define 

£(i,s):= Yl KX x ,(t) = x) 
x'eA -{ Po } 

and 

C(i,x):= Y HY x ,(t) = x), 
x 'eA' t 

where ^ := (A t -A M Po }. The processes {(^(t), g t w ) : t > 0}, {(0 Po , J*), g? '* ) = 
* > 0}, {(' w (t) : t > 0}, and {£(£) : t > 0} are then coupled. 

Let V be the first time that some particle with a label smaller than po is at the 
right-most visited site while no particle with label larger than or equal to po is at 
the right-most visited site 

V := inf [t > : £(t, g t ) > 0, Q Pom (jt, qt) = o} . 

Lemma 2. Lei A C Z, to := ({y^ : x G 74o},go) G § an d Po = sup{x : 
x G Aq}. Assume that y po = q$. Consider the corresponding coupled process 

K(Ut),qr)d<5 P0 ,Jt),q^),ti(t)y,t> 0}. Then 

q^" w =q?, t<V, 

P t PO ' qo =p?, t<V 

and 

Cs PO , qo (t) = C(t), t<V. 

Proof. It is enough to observe that before time V', none of the particles with labels 
smaller thanpo affect the dynamics of the process {({Y x (t) : x G A t }, qt) ■ t> 0}. □ 
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2.3. Auxiliary process. Let us now define a process which will be helpful to obtain 
estimates for the law of some stopping times used to define the regeneration times. 
Let n > 2 be a fixed natural number. Now let {Z x : x > 0} be a set of independent 
continuous times simple symmetric random walks of rate 2/n, such that Z x (0) = x. 
Define v\ as the first time the random walk Zq hits the site 1, 

V! := m£{t > : Z (t) = 1}. 

For m such that 2 < m < n, define v m as the first time that any of the random 
walks Zq, Z\, . . . , Z m —i hits site m, 

v m := inf{t > : sup Zi(t) = m}. 

0<i<m-l 



And for m > n, define 



v m := inf{t > : sup Zi(t) = m}. 

m—n<i<m—l 



Now define for t > 0, 



qt := m for V{ < t < 



m m+1 
j=0 i=0 



with the convention that vq = 0. In the sequel we will call {qt '■ t > 0}, the auxiliary 
process with n particles. 



Lemma 3. Whenever n > 3, there is an a > such that 



liminf > a a.s. 

t-^oo t 

Proof. For each natural n, we can construct a process having the same law as {gf '° : 
t > 0}: the last n activated particles have priority over the ones activated previously; 
nevertheless, if at a given time there are m particles from this group (of n particles), 
each one jumps at a rate 2/m. We can couple this construction with the auxiliary 
process with n particles {qt '■ t > 0} in such a way that q^°'° > qt (for the details of 
such a coupling within a similar context see [2]). Now, it is easy to check that 



lim — 

t^co t 



--: a 



a.s. 



□ 
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2.4. Regeneration times. Let us consider an initial condition ({y x : x £ Ao},qo) 
such that y Po = qo. Let also a\ and a 2 be such that < 2a\ < 02 < a- Define 

T :=M{t>0: Z Khqo (t)=q -l}, 
which is the first time that the leading particle at qo jumps backwards. Now let 

q 5 t PO ' qo =Z po , qo (t) forO<t<T 
and q S t P0 ' q ° = q S t P0 ' q ° for / > T. Define 

U := inf {t > : q^ - q < [a 2 t\ or Cs PO , qo (t, 0) = 0}, 
V := mf{t > : £(t, [aitj + g ) > 0} 

and 

D := mm{U,V}. 

The definition of U in terms of q^ P0 ' q ° instead of q^ P0 ' q ° is necessary because we 
have to avoid the possibility that before time D some particle in the configuration 
£ branches before the front q t P0 ' q ° has moved. Furthermore, the fact that 2a\ < a 2 
guarantees that for times t > \ ja\ the function [a 2 t\ is never equal to [or^J- Note 
that U, V and D are stopping times with respect to the natural filtration {J-f : t > 0} 
of the labeled zero-range reactive process. 

Let us also define the first time U and V happen after time s > 0, 

U o S := inf{/ > : qf s ><" - q s < [a 2 t\}, 

Vo9 s := inf {t > : ^ (t, [ ai t] + g 8 ) > 0}} , 

and D o9 s := min{U o 9 S ,V o S }, where w s := ({Y x (s) : x € ^4 S }, g s ). 
For each y € N, define the ^-stopping time 

N y :=mf{t>0:p t = Po + y}. 

We now define sequences of stopping times {Sk : A; > 0} and : A; > 1} as follows. 
First let S := and R := 0. Then define for k > 0, 

S'fe+i := N Rk+ i, D k+ i := D o 9 Sk+1 + Sk+i, Rk+i ■= PD k+1 , 
where we adopt the convention that pn k = 00 and S^+i = 00 in the event = 00. 
We similarly define Uk ■= U o 0s k + Sk and Vk := V o Q Sk + Sk for k > 1 . Let now 

K := ini{k > 1 : Sk < 00, = 00}, 
and define the regeneration time, 



k := S K - 
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As in [2], k is not a stopping time. Define Q, the information up to time k, as 
the completion of the cr-algebra generated by events of the form {t < k} D A, with 

The following proposition will be proved in section HI We call an element w € S 
nontrivial, if it corresponds to a configuration with at least one particle behind the 
front or at it. 

Proposition 1. For every non-trivial initial condition w(!§, 

(4) k < oo, ¥ w — a.s. 
Furthermore, 

(5) E5 00 [k 2 |C/ = oo] < oo and "E,^ [r 2 K \U = oo] < oo. 
Proposition 2. Lei F be a Borel subset of D([0, oo); $7). Then, for every nontrivial 

F w [{q K+ .-q K ,p K+ .-p K ,T- qK) - PK C K ( K + -)) e F\G] =F Sofi [{q.,p.,a-)) & F\U = ov]. 
Proof. It is enough to prove that for every B £ Q, 

F w [B,(q K+ . -q K ,p K+ . - p K , r_g K ,- Pk C k ( k + ')) e F ] 
= P[S]P 50i0 [( g .,p.,C(-))GF|i7 = O o]. 

As in [2], this can be done using lemma [2] and observing that on the event on the 
event Sk < oo, q Sk = x and p Sk = y, we have that 

whenever Uf. = V/. = oo. □ 
We can now define a sequence «i < «2 < • • • , with K\ := n while for n > 1, 

K n+ i := K n + k(u> k „ + •), 

where k(w Ku + •) is the regeneration time starting from w Kn and we set K n +i = oo on 
the event K n = oo. We call k\ the first regeneration time and K n the n-th regeneration 
time. Now define for each n > 1, the cr-algebra <5 n as the completion with respect to P 
of the smallest cr-algebra containing all sets of the form {k% < ii}(T • -n{K n < t n }f)A. 
A E Tt n ■ As in [2] , we then have the following generalization of proposition [2J 

Proposition 3. Let F be a Borel subset of D([0, oo); $7). Then, for every nontrivial 
w € S anci n > 1, 

= F 5ofi [(q.,p.,C(-))eF\U = oo]. 
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Proposit ion 4. Let w <G S. (i) Uftdcv IP^; rei, K2 — «i, K3 — «2) • • ■ ar e independent, 
and K2 — Hi, «3 — «2, • • • are identically distributed with law identical to that of n\ 
under T>s 0i0 [-\U = oo]. fraj Under F w , g.A/ci,9(« 1 +0AK2 ~ ?(k 2 +-)ak3 are 
independent, and g( Kl+ .) AK2 — g Kl , q( K2+ .) Aft3 — q K2 , . . . are identically distributed with 
law identical to that of q Kl under Fs 00 [-\U = oo] . (Hi) UnderP w , P-ak 1 iP( Ki +-)ak 2 ~ 
P/ti.P(K2+-)AK 3 ~P«2' • • • are independent, and p {Kl+ . )AK2 -p Kl ,P( K2+ .) A K 3 ~Pk 2 , ■■■ are 
identically distributed with law identical to that of p Kl under ¥s 0Q [-\U = oo] . 



3. Proof of theorems 1 and 2 

3.1. Proof of theorem 1. The proof of parts (i) and (u) of theorem 1 follow now 
using standard arguments (see for example [2] or in the context of Random Walks 
in Random Environment |12j). Indeed, using Propositions Q] and |4] we first prove 
that, a.s. 



q Kn E w [q Kl \U = oo] p Kn E w \p R \ U = oo] 
hm = =: V\ and hm = — — — L — =: V2- 

n— »oo K n K w [Ki\U = OO] n^oo K n K w [Ki\U = Oo] 

An interpolation argument then proves that a.s. 

lim — = v\ and lim — = V2- 

t— >oo t t^oo t 

Since r t = qt + Pt, this implies that a.s. 

lim = vi + v 2 , 

t— >oo t 

proving part (£) of theorem [TJ We then define Pj := p Kj+1 — p Kj — {ftj+i — Kj)v2 and 
Rj ■= P Kj+1 + <lK j+1 - Pk 3 - - - + v 2 ), and show that 

m m 

S m := ^2 P 'j and S ™ := ^2 R 3> 

3=1 3=1 

converge in law to Brownian motions with variances 



E w [(p Kl - v 2 ni) 2 \U = oo] E w [(r Kl - (vi + v 2 )ki) 2 \U = oo] 

°1 := ^ — F T77 5 an d °2 : ~ 



E^[ki|C/ = oo] K w [ki\U = oo] 

respectively. An interpolation argument can then be used to obtain the full limits 
proving that 



(6) JPfe-it-e-Hvt) and e 1/2 (r e -i t - e^i^i + v 2 )), i > 0, 

converge in law to Brownian motions with variances a\ and o~2 respectively. In 
section [5] we will prove that these variances are positive. 
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3.2. Proof of theorem 2. Let \xt be the law at time t of the exclusion reactive 
process seen from the front T- n r}t € Qq := {0, This is a Markov process with 
infinitesimal generator 

Lof(v) = w(0)(/(r_iT7 + <5 ) - f(rj)) + (1 - p)r,(0)(f(T^r,) - f(rj)) 

+ ^, v <p,\x-y\=i - v(y))(f(^ yV ) - f(r,)). 

For a local function / on fioj define 1(f) as the smallest integer I such that f(rj) does 
not depend on rj(x) if x < — I. Now define the probability measure /Uqo on Qq by the 
formula 

_ E g o[/ f y +1 /(r- ra r/ B )d S |^ = oo] 
Q °° E 5o [ki|C/ = oo] 

for TV such that N(a' — a") > 1(f)- This defines a consistent family of probability 
measures on cylinders. 

Theorem 3. hmt_»oo fit = Atoo weakly and /Uoo is invariant for the generator Lq. 

Proof. First note that for / local and w an arbitrary initial condition, 
lim^oo K w [kn+i > f( T -r t v(t))] = 0- Therefore in the decomposition 



fdjl t = E w [k N+ i < t, f(T- n f](t))] + E w [k N+1 > t, f(T- n T](t))], 

it is enough to examine the first term. Now, as in [2], by Proposition [3] 

(7) E w [k n+1 <t,f(T- rt v{t))] = I M t (du)F f (u), 

Jo 

where 

M t ([0,u]) G [t-u,t]) 



k>l 



and 



F f (u) := E$ 00 [k n <U< K N+1 , f(T-r t T)(u))\U = oo]. 

Also, as in [3], we can check the spread out assumption of the renewal theorem 
(theorem 6.2 of p3]) to show that 

ii 

limM([0,«]) 



t^oo ' E5 00 [ki|?7 = oo] 

uniformly on compacts. Finally, since \Ff\ < ||/||oo|-^i| and J F\du < oo we conclude 



that 



lim / fdfit = / fdfioo. 



□ 
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4. Estimates for the regeneration times 
Let us first state the following estimate for the stopping time U. 

Lemma 4. For every p > 1 there exists a constant C = C(p) such that for every 
initial condition w 

w w [t < u < oo] < cr p . 

Proof. The proof is based on a comparison with the auxiliary process defined in 
subsection 12.31 and large deviation estimates as in the proof of lemma 7 of [2] . □ 

Let us define V as the first time that independent continuous time simple sym- 
metric random walks created at rate 1 at the origin, reach the site \_at\. We have 
the following estimate for the stopping time V. 

Lemma 5. There is a constant C > such that for every initial condition w 

F w [t<V <oo]< exp {-Ct} . 

Proof. Let us note using the labeled construction of the zero-range reactive process, 
that we can bound the probability of event {t < V < oo} by the probability of the 
event {t < V < oo}. It is not difficult to check that there is a constant C > such 
that, P[t < V < oo] < exp {-Ct} (see [2J). □ 

From Lemmas [H and [5] we can now directly prove (see [2]), 

Lemma 6. For every p > 1, there is a constant C = C{p) such that for every initial 
condition w 

¥ w [t < D < oo] < Ct~ p . 

We continue with two important lemmas. The first one is proved using the aux- 
iliary process as in [2]. 

Lemma 7. There is a 5\ > such that for every initial condition w 

F W [U = oo]> S v 

The second lemma which follows can be proved using the inequality PiufV = oo] > 
F w [V = oo]. 

Lemma 8. There is a 82 > such that for every initial condition w 

f w [V = oo}>5 2 . 
We will also need the following lemma. 
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Lemma 9. There is a constant C > such that for every initial condition w with 
p = q = and M > I, 

^w[qt > Aft] < Cexp {-Ct} and F w [p t > Mt] < Cexp {-Ct} . 

Proof. It is enough to note that the rate at which the process {qt : t > 0} increases 
is always bounded by 1. Similarly for the process {pt ■ t > 0}. □ 

Let us now prove of Proposition [TJ Note that for every k > 1, ¥ w [k = oo] < 
^w[Dk < oo]. Taking the limit when k —* oo we obtain Part ([5]) of Proposition 
[T]is a consequence of the following lemma. 

Lemma 10. There is a constant C > such that for every p > 1 and t >0, 

p Sq , [k >t\u = oo] < cr v . 

Proof. Note that 

oo 

¥ Soo [k > t\U = oo] = ^P 5o , [5 fc >t,K = k\U = oo}. 

k=l 

From Lemmas [7] and [8l using the strong Markov property we deduce that there is a 
5 > such that for every Z > 1, 

I 

(8) P 5o > > t\U = oo] < ^ p 5o,o[Sfc > *,-K" = fc|«7 = oo] +<5- 1 (l - (5)', 

k=l 

Let < 7 < 1 and consider the event 

A k := {r Dl - r Sl < t 7 , . . . , r rifc _ 1 - r 5fc _ 1 < t 7 }. 

On Afc we have r$ k < A;t 7 . Since = % (because the initial condition is 5o,o)> if 
U = oo then r t > qt > [at\ for all t > 0. We then have on fl {{/ = oo} that 
[aS k \ < W. Hence, for t > {W + l) /a and fe < I, 

p <5o,o[* < ^fc < oo, A k \U = oo] = 

and 

P 50j0 [t < S k < oo\U = oo] < P 50j0 [t < S k < oo, A||C7 = oo] 

(9) < CE^P^k-rs, > t 7 ,5 fe < oo], 

where we have used Lemma [71 Now, using Lemmas [9] and [61 we can show that for 
every p > 1 there exists a constant C such that 

P5o>A " ^ > t\ S k < oo] < Ct-*. 
Substituting this estimate back into (JI3) and then into (JHJ), we conclude the proof. □ 
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5. NON-DEGENERACY 

To prove the non-degeneracy of the limit of the first expression in display ([6]) we 
will prove that for some a < (3 < v 2 it is true that 

(10) n , \p K1 = 1,/r 1 < ki\u = 00] > 0. 

Note that, 

F«5o,o [P«i = 1, < «i, U = 00] > ¥ SOj0 [/T 1 < St < U, D o Sl = oo]. 
The lower bound of the above inequality can be written as 

E^otlGS -1 <Si< U)E SOi0 [Uoe Sl =oo,Vo0 Sl = oo\F Sl ]]- 

Now 

E <5o,o ^ = 00, V o 6>5j = oolJFsJ = 

where in the inequality we have used Lemmas [7] and [H Thus, 

Vs ,o\P~ 1 <S 1 <U,Doe Sl = ™}> 5 1 S 2 F 5oo [p- 1 <S 1 <U}. 
Now, the probability of the right-hand term can be lower bounded by the probability 
that the particle initially at site with label 0, performs a first jump at some 
time t towards the right and then at time t + t' it branches at site 1 such that 
p- 1 < t + t' < a- 1 . Clearly this probability is positive. 

To prove the non-degeneracy of the second expression in display ([6]) this time we 
will show that for some a < (3 < v\ + v 2 it is true that 

(11) P*o>ki =hP~ 1 <ki\U = oo] >0. 

But note that {r Kl = 1} = {p Kl = 1}. Therefore inequality (fTTj) follows from (fTUj) . 

6. AN ALTERNATIVE DEFINITION OF REGENERATION TIMES 

6.1. Asymmetric simple random walk. Let us show through a simple example 
how to prove the independence of the increments of the regeneration times, with- 
out defining them through a sequence of stopping times. What we present here is 
essentially contained in the paper of Kesten [6]. 

Consider a discrete time asymmetric simple random walk {X n : n > 0}, which at 
each site jumps to the right with probability q and to the left with probability p, 
with q > p. Define the random time 

K := min{n > : minXfc > maxX,-}, 

k>n j<n 

with the convention that maxj<o Xj = — 1. In words, k is the first time the random 
walk visits a new site, without never afterward moving to the left of such a site. 
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Let us define the cr-algebra T K of k as the one generated by events of the form 
{k = n} n {X\ = xi, . . . , X n = x n } for n > and x\, . . . ,x n in 7L d . If we call P the 
law of the random walk {X n }, starting from 0, the following equality is satisfied for 
every subset A G Z N : 

Lemma 11. 

(12) P [X K+ . - X K G A\T K ] = P [X. G A\k = 0] . 

Let us prove this equality using directly the Markov property. For convenience, 
we define Y n := s\xp k<n Xk, for n > 0. Note first that 



(13) {k = n} = { inf X m > X n } n { sup X p = X n - 1} n B n , 

m >n l<p<n-l 

where 



Bn '■= {VA; < n : there exists a such that k < j < n and < 1^ 
or there exists a j, such that < j < k and Yj = X^} . 

Let us remark that B n € J- n , where {J- n : n > 0} is the natural filtration of 
the random walk {X n : n > 0}. Furthermore, the condition q > p, ensures that 
P[« = n] > 0. Defining C n := {sup 1<p<n „ 1 X p = X n — 1} PI B n G we now have 

P [X K+ . - X K e A|k = n, Xi = xx, . . . ,X n = x n ] 
_ P [X K +. — X K g A, inf m > n X m > X n , C n ,X\ = x±, , , , , X n = x n ] 
P [inf m>n X m > X n , C n ,Xi = x\ 1 . . . , X n = x n ] 

On the other hand, by the fact that C n G J- n , the Markov property, and translation 
invar iance, we have that 



P 



E 



X K+ . X K G j4, inf X m > X n , C n , Xi — xi, . . . , X, 

m>n 



n 



L{C n ,X 1 =X 1 ,...,Xn=X n } 



P 



X n +. — X n G A, inf X m > X n 



m>n 



(14) = P[C n ,Xi =xi,...,X n = x n ]P[X. G A, inf X m > 0]. 

m>0 



Choosing A as the whole space in the previous development, we conclude that 



P 



inf X m > X n , C n , X\ — x\, . . . , X n — x n 

m>n 



P[C n ,X x = Xl ,...,X n = x n ]P[ inf X m >0], 

m>n 



which combined with (|14p completes the proof of (jl2[) . 
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6.2. Exclusion reactive process. Let us now show how the approach presented 
in the context of a symmetric simple random walk can be implemented for the 
exclusion reactive process to define a version of the regeneration times. To define 
the regeneration times we will consider the holes (empty sites) of the process as 
second class particles. 

We will construct the exclusion reactive process associating to each bond con- 
necting two nearest neighbor sites independent Poisson processes each one of rate 1. 
Let us assume that the initial condition is of the form (0,7/(0)), where 7/(0) is any 
nontrivial configuration of particles. If at a given time the rightmost visited site is 
r and x + 1 < r, whenever the Poisson clock connecting sites x and x + 1 rings, the 
state at sites x and x + 1 is interchanged: if the state of the process was rj, this is 
changed to a X:X+ irj. In this case the front stays at r. If at a given time the rightmost 
visited site r is occupied and the Poisson clock connecting sites r and r + 1 rings, 
with probability p a particle is created at site r + 1 and with probability 1 — p the 
particle at r jumps to site r + 1. In both cases the front advances one step to r + 1. 
Define t\ as the first time the front advances one step creating a new particle by 

T\ := inf{t > : r t > and r\{r t — 1, t) = r](r t , t) = 1}. 

In general for n > 2 define recursively T n as the first time after time r n _i that the 
front advances one step creating a new particle by 

T n : = inf{t > r n _i : r t > r Tn _ x and rj(r t - l,t) = v( r t,t) = 1}. 

For each n > 1, at the time T n we will fill up the holes defined by the configuration 
7/(r n ) with particles: the original particles in this configuration are first class with 
respect to the new ones, which we will call holes. After time r„ holes can never 
activate a new particle. 

Let us now adopt the convention that whenever the Poisson clock of a bond 
corresponding to two occupied sites rings, the corresponding particles or holes are 
interchanged. We can then label each particle and each hole, following its trajectory. 
We will assign the label to the particles initially at r or to the left of r. We assign 
the label 1 to the particle activated at time t\. In general for n > 2, we assign the 
label n to the particle activated at time r n . Similarly, we assign the label to all the 
holes initially at r or to the left of r, and to the new holes created before time t\. 
We assign the label 1 to the holes created after time t\ but before time T2, and for 
n > 2 we assign the label n to the holes created after time r n but before time T n+ \. 
Let us call T the state space of this process consisting of the ordered pairs (r,rj), 
where r is an integer representing the position of the front and rj is a configuration of 
labeled particles and holes at sites . . . , r — 1, r. Let us call Q the law of this process 
in the corresponding Skorohod space. We now define for each n > 1, the stopping 
time D n as the first time after time r n that some of the particles or holes with labels 
strictly smaller than n (thus, excluding the foremost particle, created at time r„) is 
at the front 



D n := inf{t > : t/(t„ + 1, r Tn+t ) = 1}. 
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Here, at time t > 0, f]{r n + 1) is the particle-hole count of those particles or holes in 
ij( T n + 1) with labels strictly smaller than n. Now define the first regeneration time 

K\ := inf{r m : m > l,D m = oo}. 

We then define recursively for n > 2, 

K n := inf{r m : r m > K n -i,D m = oo}. 

Let us call 5q G T any initial condition with rightmost visited site r = 0, one 
particle with label at and none elsewhere, and one hole at each site x < 
each one with label —1. Given this initial condition, we define Do as the first 
time that one of the holes with label —1 is at the front. Call Q\ the information 
up to time K\. We have the following proposition corresponding to Proposition 
El We will in general call rj the particle count corresponding to a state 77 and 
W := {(r, ij):r£Z,i|e{0, 

Proposition 5. Let F be a Borel subset of D([0, 00); Q'). Then, for every nontrivial 
rj G T, which has only particles with label and holes with label 0, 

Q v [(r K1+ . - r Kl ,T_ rKiV '( Kl + •)) G F\g x ] = Qs [(r.,v(-)) G ^lA) = 00], 

where 7/ is £/ie particle count of those particles created at time k\ or after. 

Proof. Note that the events A n {k\ = r n } with A G T Tn for some n > 1, generate 
the ex-algebra £1. Therefore, it is enough to prove that 

(15) Q v [(r Kl+ . - r K1 ,T-r K1 rf(Ki + •)) G -F, A«i = r n ] 

= Q<5 [( r -' r ?(-)) G F|L> = oo]Qs [A,Ki = r n ], 

for each natural n and ^4 G •?>„• Now, note that {ki = r n } = {L> n = 00} n {D\ < 
00} . . . n {-D n -i < 00}. But note that for every 1 < j < n, {D n = 00} n {Dj < 
00} = {D n = 00} n {Dj < r n }. Therefore 

{Kl = T„} = {.Dn = OO} PI {£>! < T n } . . . PI {D n -l < T n }. 

Defining B := {D\ < r n } . . . n {D n -\ < r„}, it follows that the right-hand side of 
(fT5l) equals 

Qvl( r Ki+- ~ r «D r -r K1 '? / («l + •)) G = T n ] 

= Qvi( r r n + - ~ r r n ,r-r Tn V'(Tn + ")) G F, A, 5, D n = Oo] 

= Q v [A,B,Q v [(r Tn+ . - r Tn ,T- rTn i(T n + ■)) eF,D n = oo[.7>J]. 

Now, by translation invariance and the strong Markov property, Qr)[( r T n +- — 
r Tn ,T-r Tn i(T n + •)) G F,D n = co\T Tn } = Qs [(r.,Tj(-)) G F, D = oo]. Hence, 



(16) 



Qvl( r Ki+- - r Ki,T-r K1 v'(^l + •)) G F, A, Kl = T n ] 

= Q V [A,D 1 <r n ,...,D n _i <T n ]Q So [(r.,ri(-)) €F,D = 00]. 
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Choosing F = Q' in the above equality we see that, Q V [A, k\ = r n ] = Q V [A, D\ < 
r n , . . . , -D n -l < T n ]Qs [Do = oo]. Substituting this back into (fTBT) we conclude the 
proof of the proposition. □ 
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